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LECTURE 27. COMPLEX SOLUTIONS AND THE FUNDAMENTAL MATRIX 



Complex eigenvalues. We continue studying 
(27.1) y' = Ay, 

where A = {aij) is a constant n x n matrix. In this subsection, further, A is a real matrix. When A 
has a complex eigenvalue, it yields a complex solution of (27.1). The following principle of equating 
real paris then allows us to construct real solutions of (27.1) from the complex solution. 

Lemma 27.1. If y{t) = d{t) + if3{t), where a{t) and (3{t) are real vector-valued functions, is a complex 
solution of (27.1), then both d{t) and P{t) are real solutions of (27.1). 

The proof is nearly the same as that for the scalar equation, and it is omitted. 

Exercise. If a real matrix A has an eigenvalue A with an eigenvector v, then show that A also has 
an eigenvalue A with an eigenvector v. 



1 - A 1 
-4 1 - A 



has two 



Example 27.2. We continue studying A = . Recall that paW 

complex eigenvalues 1 it 2i. 

If A = 1 + 2i, then A — \I = has an eigenvector ■ The result of the above 

exercise then ensures that ( 1, . ) is an eigenvector of the eigenvalue A = 1 — 2i. 



In order to find real solutions of (27.1), we write 

(l+2^)t cos2t \ t f sm2t 

\2iJ \-2 sin 2tJ ^ \2 cos 2t 

The above lemma then asserts that e* ( q^^-^I. ) and e* { n^^'^lA are real solutions of (27.1). 

\ ^ Sill j \ ^ COS ZiXj I 

Moreover, they are linearly independent. Therefore, the general real solution of (27.1) is 

t f cos 2t \ I f sin 2t 
\^-2sin2y +^2*^ \^2cos2t 

The fundamental matrix. The linear operator Ty := y' — Ay has a natural extension from vectors 
to matrices. For example, when n = 2, let 

yiA - ffA f yiA - (91 



y2i) v/2/ ' \y22) \g2 

Then, 

J, ( yii yi2\ ^ (h ai 

\y21 ^22/ 92, 

In general, if A is an n x n matrix and Y = {yi ■ ■ ■ y„) is an n x n matrix, whose j'-th column is yj, 
then 

TY = T{yi ■■■ yn) = {Tyi ••• Tyn). 
In this sense, y' = Ay extends to Y' = AY. 
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Exercise. Show that 

T{U + V) = TU + TV, T{UC) = {TU)C, T{Uc) = {TU)c, 
where U, V are n x n matrix-valued functions, C is an n x n matrix, and c is a column vector. 

That means, T is a linear operator defined on the class of matrix-valued functions Y differen- 
tiable on an interval /. The following existence and uniqueness result is standard. 

Existence and Uniqueness result. . If A{t) and F{t) are continuous and bounded (matrix-valued 
functions) on an interval to G /, then for any matrix Yq then initial value problem 

r = A{t)Y + F{t), Y{to) = Yo 

has a unique solution on t G /. 

Working assumption. A{t), F{t), and f{t) are always continuous and bounded on an interval 
t E /. 

Definition 27.3. A fundamental matrix of TY = is a solution U{t) for which |C/(to)| / at some 
point to- 

We note that the condition \U{to)\ / implies that \U{t)\ / for all t G /. We use this fact to 
derive solution formulas. 

As an application of U (t), we obtain solution formulas for the initial value problem 

y' = A{t)y + f{t), y{to)=yo. 

Let U{t) be a fundamental matrix of Y' = A{t)Y. In the homogeneous case of f{t) = 0, let 
y{t) = U{t)c, where c is an arbitrary column vector. Then, 

y' = U'c = {A{t)U)c = A{t){Uc) = A{t)y, 

that is, y is a solution of the homogeneous system. The initial condition then determines c and 

c = U-\to)yo. 

Next, for a general f{t), we use the variation of parameters by seting y{t) = U{t)v{t), where v 
is a vector-valued function. Then, 

y' = [Uv)' = U'v + Uv' = A{t)Uv + Uv' = A{t)y + Uv'. 
Hence, Uv' = f{t) and 

y{t) = U{t) I U-\t)f{t)dt. 

Liouville's equation. We prove a theorem of Liouville, which generalizes Abel's identity for the 
Wronskian. 

Theorem 27.4 (Liouville's Theorem). IfY'{t) = A{t)Y{t) on an interval t e I, then 
(27.2) \Y{t)\' = tiA{t)\Y{t)\. 

Proof. First, if \Y{to)\ = at a point to e /, then \Y{t)\ = for all t £ I, and we are done. We 
therefore assume that \ Y{t)\ / for all t G I. 
Let Y{to) = / at a point to- That is, 

Y{to) = (yi(to) • • • y„(to)) = {Ei E2 ■ ■ ■ En). 

Here, Ej are the unit coordinate vectors in M", that is, the n-vector Ej has 1 in the j-ih position 
and zero otherwise. 
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We use the derivative formula for the determir\ant 

|y(t)r=|det(yi(t)...y„(t)) 

= det(yi(t) y2{t) ••• + det(yi(t) ••• + • • • + det(yi(t) ••• y'^{t)). 

This formula is based or\ the Laplace expansion formula for determinant, and we do not prove it 
here. Since 

y'jito) = A{to)yj{to) = A{to)Ej = Aj{to), 
where Aj{t) is the jth column of A{t), evaluating the above determinant formula at t = to we 
obtain 

\Y{to)\' = det(Ai(to) E2 ■■■ En) + det(Si ^2(^0) ■ ■ ■ En) + ■ ■ ■ + det(Ei E2 ■ ■ ■ An{to)) 

= an (to) + 022(^0) H h ann{to) = tryl(to)- 

Thus, (27.2) holds at to- 

In general, let C = Y{tQ)-^. Then U{t) = Y{t)C satisfies 

U' = A{t)U, [/(to) = /. 

Therefore, by the argument above |f7(to)|' = tr^(to)|C/(to)| = tryl(to). Since 

|(iy(t)ci) = j^mmc\) = \Ymc\, at t = to, 

it follows that tr^(to) = |y(to)ny(to)|^^. Since to is arbitrary, the proof is complete. □ 
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